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Abstract. Spectroscopic radial velocities for several nearby open clusters suggest that spectra of (especially earlier- 
type) rapidly rotating stars are systematically blueshifted by 3 km s~ or more, relative to the spectra of slowly 
rotating ones. Comparisons with astrometrically determined radial motions in the Hyades suggests this to be an 
absolute blueshift, relative to wavelengths naively expected from stellar radial motion and gravitational redshift. 
Analogous trends are seen also in most other clusters studied (Pleiades, Coma Berenices, Praesepe, a Persei, IC 
2391, NGC 6475, IC 4665, NGC 1976 and NGC 2516). Possible mechanisms are discussed, including photospheric 
convection, stellar pulsation, meridional circulation, and shock-wave propagation, as well as effects caused by 
template mismatch in determining wavelength displacements. For early-type stars, a plausible mechanism is 
shock-wave propagation upward through the photospheric line-forming regions. Such wavelength shifts thus permit 
studies of certain types of stellar atmospheric dynamics and - irrespective of their cause - may influence deduced 
open-cluster membership (when selected from common velocity) and deduced cluster dynamics (some types of 
stars might show fortuitous velocity patterns). 
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1. Introduction 

Wavelength displacements in stellar spectra are affected 
not only by stellar radial motion causing a Doppler shift, 
but also by line asymmetries and shifts caused by, e.g., 
atmospheric pulsation, surface convection, stellar winds, 
and gravitational potential. However, for solar-type stars, 
the net effect of convective blueshifts and gravitational 
redshifts, as predicted from hydrodynamic model atmo- 
spheres, does not exceed 1 km s _1 (e.g., Dravins 1999). 

For a few nearby open clusters, it has been possible to 
determine accurate radial motions from astrometry only, 
not invoking spectroscopy. Comparing such astrometric 
radial velocities to spectroscopic data for the Hyades, in- 
dicates not only the existence of such convective lineshifts, 
but also an increasing blueshift of stellar spectra, reaching 
3-5 km s^ 1 for rapidly rotating stars earlier than about 



F5 V: Madsen et al. (20021) and Dravins (2003). Since 
there do not yet exist any detailed hydrodynamic models 
of such earlier-type stellar atmospheres from which syn- 
thetic spectral-line shifts could be reliably computed, it is 
not yet possible to theoretically verify the reality of these 
suggested trends. 

In order to examine this effect, we have analyzed spec- 
troscopic data for as large a dataset as possible, using 



the open-cluster database WEBDA 1 . Although its data 
sets are somewhat inhomogeneous, this permits a study 
of relative wavelength displacements for ten clusters, as 
function of stellar spectral type and of rotational velocity. 

Five of these clusters have astrometric radial veloci- 
ties deduced from Hipparcos data (Madsen et al. I2()()2|l : 
Hyades, Pleiades, Coma Berenices, Praesepe and a Persei, 
but only for the Hyades (and perhaps also Coma) is the 
accuracy adequate to permit studies of not only relative 
but also absolute wavelength shifts. For stars not mea- 
sured from Hipparcos, their astrometric radial velocities 
are estimated by assuming they follow the cluster's com- 
mon space motion. 

Spectroscopic radial velocities for the other five clus- 
ters (IC 2391, NGC 6475, IC 4665, NGC 1976 and NGC 
2516) permit studies of relative wavelength shifts relative 
to their mean. This is possible for clusters of small angular 
size whose stars share the same velocity vector, so that the 
variation across the sky of the radial motions is negligible 
compared to the shifts one is searching for. 

All ten clusters have metallicities close to solar. 
Their ages range from 13 Myr (NGC 1976) to 660 Myr 
(Praesepe). 
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2. The 10 clusters 

Fig. d shows absolute wavelength shifts of stellar spectra 
in the Hyades, plotted as spectroscopic minus astrometric 
radial velocities. A slight trend of a somewhat increased 
blucshift (by about 1 km s _1 ) when going from K-type to 
mid-F type stars (B — V from about 1.1 to 0.5) is expected 
due to the more vigorous convection in F stars, causing 
increased convective blueshift. However, this trend contin- 
ues and becomes enhanced for even hotter stars, a region 
where not many detailed hydrodynamic models have yet 
been developed. Many of these stars are also rapidly rotat- 
ing, making it difficult to disentangle whether the depen- 
dence is upon spectral type or/and on rotation, or perhaps 
even on something else. A similar dependence on rotation 
is suggested also from data for some cool Hyades M-stars: 
Fig. The trend seen in the Hyades is also present in 
Praesepe and possibly in Coma Berenices. In a Persei the 
same is seen except for the most rapid rotators (all of 
spectral type B). For the Pleiades, the trend for earlier 
spectral types is at best weak, while its M-stars do show 
a clear correlation, noted also by Terndrup et al. (2000); 
see their Fig. 3. The error bars in Fig. include both the 
stated spectroscopic errors and the part of the astrometric 
error coming from the internal velocity dispersion in each 
cluster. For detailed discussions on the astrometric error 
budget, see Lindegren et al. H2000[) and Madsen (2003). 

The rotational broadening of spectral line profiles 
makes the spectroscopic determination of wavelength shifts 
difficult in rapidly rotating stars, and even the errors of 
several km s _1 quoted in the literature may be optimistic. 
Therefore, one should not draw conclusions from individ- 
ual data points for rapid rotators in, e.g., Fig. 1, but only 
from their averages and from the general trends. The sta- 
tistical scatter may well account for some stars showing 
apparent redshifts. 

The trend of increasing blueshift with higher V sin i 
can be seen in IC 2391 and NGC 2516, partly in IC 4665, 
maybe in NGC 1976, but not in NGC 6475 (Fig. UJ. In 
contrast with a Persei and maybe the Pleiades, NGC 1976 
and IC 4665, the most blueshifted stars in IC 2391 and 
NGC 2516 happen to be of B-type. It must be noted that 
all these five clusters generally lack mid- and late A stars, 
for which the trend is well visible in the Hyades. There is 
no apparent correlation between the suggested trends and 
cluster age. For instance, IC 2391 and IC 4665 are coeval 
(45 Myr), while the Pleiades (130 Myr) and NGC 2516 
(110 Myr) have about the same age, but show different 
trends. 

The compiled data sets of spectroscopic radial veloci- 
ties and projected rotational velocities V sin i are quite in- 
homogeneous, especially for the early spectral types. Since 
the number of sources is huge, we only give some refer- 
ences to observations or compilations for the interesting 
part of Figs. 1 and 2 (i.e. where V smi shows a trend to- 
wards blueshifted spectra). The Hyades set is a mixture of 
the compilations by Perryman et al. I|1998i and Glebocki 
& Stawikowski (2000), and the observations by Griffin et 



> 

i 



-5 - 




Gravitational redshift 

V sin i o 10 km s"' 
Vsin i Q 100 km s" 1 



0.0 



0.5 



1.0 
B- V 



1.5 



2.0 



Fig. 1. Wavelength displacements of stellar spectra in 
the Hyades. Differences between spectroscopically mea- 
sured radial velocities, and astrometrically determined 
radial motion are plotted as function of both B — V, 
and projected rotational velocity Tasini (proportional 
to symbol area). Spectra of rapidly rotating early- type 
stars appear systematically more blueshifted than those 
of slower-rotating later-type ones. Only points with nom- 
inal differential velocity errors < 3 km s _1 are plotted. 
In the absence of effects in stellar atmospheres or due 
to spectroscopic techniques, the naively expected spectral 
displacements should equal the marked line which esti- 
mates the gravitational redshift for main-sequence stars. 
Corresponding spectral types for main-sequence stars are 
marked at the top. 



al. ljHjgg|) and Reid k Mahoney l(2"UUU)l . The Pleiades set 
mainly comes from Liu et al. (TMT)) . Morse et al. (1991), 
Raboud k Mermilliod fiMfy . Queloz et al. ifTMSjl and 

Terndrup et al. I|2000|l . Coma Berenices data are largely 
from Abt k Willmarth (1999). Praesepe data are chiefly 

mm 



from McGee et al. (Q267), Rachford : 199,§ and Abt & 
Willmarth (1999). Data for a Persei primarily come from 
Kraft fTT^l . P etrie k Heard (JMOI), Mill ward fc Walker 
IPSSI) . Prosser HTWIk and Morrell & Abt (1992J). For IC 

2391, the data are mainly from Levato & Garcia {1984) 

mm 



and Levato et al. ( 1988J); for NGC 6475 mainly from Abt k 
Jewsbury QHRty . Abt et al. ifiTTTnjl and Gieseking fT577)l . 
and for IC 4665 predominantly from Abt k Snowden 
Abt & Chaffee ifTMT) and Morrell k Abt l(TMT)l . 



NGC 1976 data are for the most part from Abt et al. 
JMOI), Smith et al. fTOgty . Morrell k Levato ijUMj) and 
Abt et al. (JTMTJ), while those for NGC 2516 mainly orig- 
inate from Abt et al. fiMfy . Abt k Levy 1(1572)1 . Herwig 
tlUHTJl and Gonzalez k Lapasset l|2000|) . 
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Fig. 2. Differences between spectroscopic and astrometric radiai vefocities versus steilar rotationai velocity Vsmi 
(top), and relative spectroscopic velocities for clusters without such astrometric data (bottom). For each cluster, the 
values are normalized to have AV = correspond to the median radial velocity of stars with Vsmi < 60 km s _1 . The 
Hyades: black points are from Madsen et al. (2002); open circles are M dwarfs (known binaries excluded) from Reid & 
Mahoney (2000); asterisks are other assumed cluster members. The Pleiades: Open circles are M dwarfs from Terndrup 
et al. (20 0"U|l . For the Hyades, stars with errors > 3 km s _1 are omitted; for other clusters those with > 10 km s — 1 . 
Those with no errors in the literature were assigned the value 2 km s _1 . 



Instead of open clusters, binaries may be used to study 
lineshifts in stars (e.g. Pourbaix et al. I2002fl . However, 
studies of visual-spectroscopic binaries (Pourbaix [5000) , 
including such with one component being a rapidly rotat- 
ing A star, have not revealed trends of blueshifts correlated 
with stellar rotation, but the effect could be hidden in the 
errors of the estimated binary parameters. 

3. Spectrum shifts due to atmospheric dynamics 

Although data for different clusters are thus somewhat 
scattered, there still remains the suggestion of a common 
trend, such that more rapid stellar rotation (and early 
spectral type) correlates with increased spectral blueshift. 
In this Section, we now look at some plausible mechanisms 
that might explain such shifts. 

3.1. Convection 

Convective blueshifts of spectral lines occur in solar-type 
stars due to the correlation between upward motions of 
photospheric gas, and their local brightness. Rising and 
blueshifted elements of hot (bright) gas contribute a larger 
number of photons than the same gas when it has cooled 
off (is darker) and is sinking (thus redshifted). The effect 
for the Sun is typically 300-400 m s -1 , increasing to about 



1 km s 1 in hotter F-type stars (Allende Prieto et al. 
I2()()2bl) . 

The effects suggested by Figs. ^ and |21 however, are 
of greater magnitude, and are unlikely to be explained by 
convective blueshifts. This assessment is based on a the- 
oretical grid of 55 two-dimensional hydrodynamic model 
atmospheres (Ludwig et al. I1999[) for A- to K-type main- 
sequence stars with solar metallicity. Fig.[3]shows the run 
of the horizontal and vertical component of the convec- 
tive velocity at Rosseland optical depth unity, as function 
of B — V . For the conversion from effective temperature 
to color, a semi-empirical relation in Alonso et al. l|1996f) 
was used. Some scatter is apparent: the large velocities in 
six models around B — V = 0.3 are due to radial oscil- 
lations which are not of convective origin, and maximum 
convective velocities of about 2.2 times solar are found at 
B — V = 0.33 (7000 K). For a given stellar surface contrast, 
the velocity scaling should be indicative of also the scaling 
for line shifts. For mid F-type dwarfs some enhancement 
occurs due to an increased granular temperature contrast 
but even then the convective blueshift would not exceed 
about 1 km s , at least not in the optical wavelength 
region. Velocities and temperature contrasts quickly drop 
towards A-type objects, reducing the shifts. 

The convective blueshift so far discussed is not corre- 
lated with stellar rotation. However, a dependence could 
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arise if rapid rotation would somehow modify the convec- 
tive flows under the influence of Coriolis and centrifugal 
forces, enhancing the blueshift. 

The horizontal scale of granulation seems to be pro- 
portional to the pressure scale height at optical depth 
unity (Freytag et al. H997[) . The velocities of Fig. lead to 
convective turn-over time scales for A- to G-type dwarfs 
within a factor 2 of the solar value (« lOmin). Even 
at rotation rates close to break-up, turn-over time scales 
remain more than one order of magnitude shorter than 
the rotational period, rendering Coriolis forces rather in- 
effective for modulating the convective surface structures. 
Centrifugal forces reduce the effective surface gravity, pro- 
ducing more "giant-like" conditions. At given effective 
temperature, lower-gravity models tend to produce larger 
shifts, but observed velocity spans of line bisectors in gi- 
ants hardly exceed 2 km s _1 (Gray & Nagel 1989, Allende 
Prieto et al. I2()()2a|l . Since large lineshifts can be expected 
to be accompanied by also large line asymmetries, it ap- 



pears improbable that lineshifts of 3-5 km s 
produced by effects of rotation on convection. 
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Fig. 3. Convective velocities from a grid of 55 main- 
sequence hydrodynamical model atmospheres (Ludwig 
et al. 1199 9). Shown is the root-mean-square horizontal 
(open) and vertical (solid circles) convective velocity at 
Rosseland optical depth unity, as function of B — V . The 
corresponding effective-temperature scale is at the top. 



3.2. Pulsation 

Gas motions where upward velocity correlates with in- 
creased brightness may cause blueshifted spectra, irre- 
spective of the scale of the gas flows. Quite a few among 
the rapidly rotating earlier-type stars showing spectral 
blueshifts reside in the area of the Hertzsprung-Russell 
diagram populated by Delta-Scuti-type and similar vari- 
ables. The sign of the effect appears correct in the sense 
that maximum outward motion in these pulsating stars oc- 
curs near the phase of stellar maximum brightness; thus 
an average over a pulsation cycle should lead to somewhat 
blueshifted lines. 

Although the number of [known] pulsating Delta Scuti 
stars in open clusters is not especially great (Antoncllo 
& Pasinetti Fracassini 119981 Rodriguez & Breger l2T)TTH) . 
several are rapidly rotating and might influence our sta- 
tistical trends. Detailed studies of the radial-velocity and 
line-shape variations have been made for a few larger- 
amplitude Delta Scuti stars (e.g., Dravins et al. 119771 Yang 
et al. H987[) . but larger pulsational amplitudes are confined 
to the slow rotators (Breger 2000). Since typical ampli- 
tudes in Delta Scuti stars are modest (maybe 3 km s -1 ), 
and typical photometric amplitudes only some per cent, 
these appear inadequate for causing a statistical bias of 
the magnitudes suggested by Figs. ^ and |2] 

3.3. Meridional circulation 

Wavelength shifts may offer a possibility to study large- 
scale stationary flows across stellar surfaces, such as 
meridional circulation. In a rotating star, meridional flows 
are generated with - in the simplest case - rising motions 
at the poles and sinking at the equator (Sweet I195()|l . 
The classical theory predicts very slow motions, occur- 



ring on the Kelvin-Helmholtz time-scale of the star, im- 
plying slow flows even in fast rotators. However, the clas- 
sical theory becomes questionable at the flow boundaries 
(for a discussion see Tassoul 20001 due to divergent so- 
lutions towards infinite flow speeds. Hitherto, this prob- 
lem has not been clarified, not excluding the possibility 
that meridional flows in the atmosphere could be faster 
than at first expected. In our case, the measured spectral- 
line displacements originate mainly from the upper stellar 
photosphere, while the classical theory is about motions 
in stellar interiors. Issues include whether the very much 
lower densities in stellar atmospheres (as compared to the 
interiors) somehow might cause more rapid flows there (to 
perhaps conserve momentum) , whether there are counter- 
rotating "cells" near stellar surfaces (Tassoul & Tassoul 
1983; their Fig. 1), or other dynamic issues connected with 
the relatively thin line- forming layers. 

Related large-scale flows at the stellar surface are ther- 
mally driven "winds" : rotation leads to gravity darken- 
ing, i.e. a non-uniform temperature distribution, with the 
stellar poles hotter than the equatorial regions. A local- 
ized heating of an atmosphere was studied by Showman & 
Guillot (2002$ in the context of driving (zonal) winds on 
irradiated extra-solar giant planets. They find wind speeds 
up to 2 km s _1 for temperature contrasts similar to those 
expected on rapidly rotating stars. 

The radial-velocity component of a meridional circu- 
lation pattern is much smaller than its tangential compo- 
nent (e.g. Tassoul 2000). The tangential velocity is - again 
in the simplest case - directed from the poles towards the 
equator of a rotating star. Such a flow seen equator-on 
would spectroscopically produce a net blueshift, seen pole- 
on a net redshift. Stars having a large V smi are likely 
observed nearly equator-on, and for such stars one would 
expect blueshifts. However, even if one assumes little can- 
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cellation by the averaging of the velocity field over the sidcrcd for explaining the intrinsic blueshifts in rapidly ro- 

visible stellar hemisphere, the velocities needed to pro- tating early-type stars. Further observational tests should 

duce the observed lineshifts have to be comparable to include analyses of differential wavelength shifts within 

the sound speed. Such flow velocities appear difficult to the spectra of individual stars, both between their weaker 

sustain under the dissipative influence of shock formation and stronger lines, and between the flanks and wings of 

(e.g., Landau & Lifschitz I1982|l . stronger ones. 



3.4. Shock-wave propagation 

Despite the fact that full 3-dimensional hydrodynamic 
models with ensuing spectral-line synthesis have not yet 
been developed for A-type and hotter stars, various 2- 
dimensional atmospheric models are becoming available. 
One sequence of such models, at temperatures near the on- 
set of convection, are being developed by Holweger ( 2003 ) . 
Models for main-sequence A-type stars of 8200 and 9000 
K effective temperatures show that the atmosphere is 
highly dynamic. It moves up and down in a pulsation- 
like manner and is traversed by numerous shocks which 
pass up through the photosphere at supersonic velocities 
of typically 20-30 km s -1 . Rising shocks (moving towards 
the observer) appear blueshifted; they are accompanied 
by density enhancements and plausibly must contribute 
blueshifted spectral-line components. Given their great 
amplitudes, it does not seem unreasonable that a statis- 
tical bias of the blueshifted profiles could well amount to 
several km s^ 1 , the amounts we have been discussing. 

In the absence of a detailed treatment of line forma- 
tion in these models, various weighting functions can be 
applied on the vertical velocity field to get an indication of 
likely wavelength shifts for the spectral lines. By weight- 
ing the local vertical velocity field with the local internal 
energy Ei (a rough proxy for spectral-line excitation), a 
displacement is found, mimicking spectral blueshifts of 1- 
2 km s _1 (Holweger 

A further factor is that the corresponding radial- 
velocity signatures probably become better visible in 
rapidly rotating stars. The wavelength shift of any stel- 
lar spectrum is normally determined through a cross- 
correlation of the observed spectrum with some spectral 
template, and then the main weight comes from the flanks 
of the sharper stellar spectral lines. For slowly rotating 
stars, many photosphcric lines (originating in the lower 
atmosphere) contribute to the radial-velocity signal. In 
rapid rotators, however, the spectral broadening causes 
most weak lines to become very shallow and not con- 
tributing much to the radial-velocity signal which is then 
weighted towards the wavelength positions of the flanks 
of the stronger lines (such as the Balmer lines from hy- 
drogen). Those typically originate in the upper photo- 
sphere, at those greater atmospheric heights where the 
rising Shockwaves have become more fully developed, and 
where one could expect their blueshifted contributions to 
be more significant. 

Even if, obviously, any genuine comparison with theory 
must eventually include detailed line-profile synthesis, we 
find this scenario the most plausible among the ones con- 



4. Systematic shifts due to template mismatch 

Systematic wavelength shifts could also result from in- 
strumental effects in the process of measuring and de- 
ducing radial velocities. Radial-velocity instruments often 
use a cross correlation template taken from a slowly ro- 
tating cooler star. However, for a rapidly rotating hotter 
star, such a template may have a zero-point mismatch in 
wavelength since the rotational broadening blends differ- 
ent spectral-line groups than those in cooler stars (e.g., 
Morse et al. ITMT1 Verschueren et al. Griffin et al. 

Even if effects of this type are present to some degree, 
we find it unlikely that they could amount to the several 
km s -1 suggested by observations (identified effects rather 
appear to lie on levels of 0.5-1 km s _1 ). 

5. Effects on open-cluster studies 

Irrespective of their cause - physical or instrumental - 
the existence of systematic blueshifts in some fraction of 
open-cluster stars can influence precise studies of cluster 
membership and dynamics. For stars whose membership 
in a cluster is doubted, a probability of membership is 
normally computed, taking into account the observed de- 
viation of the star's motion from that of the cluster as a 
whole. Obviously, if for some groups of stars such mea- 
sures are systematically shifted due to other effects, these 
stars may erroneously be rejected as cluster members (or 
non-members included). 

Also biases may enter for any kind of studies of cluster 
dynamics that are based upon spectroscopic radial veloci- 
ties. The effects are likely to be greatest for groupings with 
many rapidly rotating stars, in particular young clusters 
and OB-associations. 

6. Conclusions 

Absolute and relative wavelength shifts permit certain 
studies of stellar atmospheric dynamics. For some nearby, 
well-studied open clusters, indications are that rapidly ro- 
tating stars (especially those of early spectral types) have 
their spectra blueshifted by 3-5 km s -1 . Having searched 
for mechanisms that might produce such shifts, we con- 
clude that one plausible candidate appears to be upward 
propagating shock waves, now seen in hydrodynamic simu- 
lations of hotter stars. However, other mechanisms cannot 
be excluded, and will probably be required to explain the 
apparently analogous trends seen also for rapidly rotat- 
ing M-type dwarfs. Irrespective of their causes, such shifts 
may influence assigned cluster membership for any star, 
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and affect the deduced cluster dynamics. The phenomenon 
may also imply that rapidly rotating field stars could have 
underestimated radial velocities. 

Future analyses of these phenomena should include, on 
the observational side, large and homogeneous data sets, 
and also separate wavelength measurements for different 
classes of spectral lines. On the theoretical side, the pri- 
mary need is for synthetic spectral line profiles computed 
from hydrodynamic model atmospheres. 
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